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Abstract 

In this paper we analyse the integrabiHty of a dynamical system describ- 
ing the rotational motion of a rigid satellite under the influence of gravitational 
and magnetic fields. In our investigations we apply an extension of the Ziglin 
theory developed by Morales-Ruiz and Ramis. We prove that for a symmetric 
satellite the system does not admit an additional real meromorphic first integral 
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except for one case when the value of the induced magnetic moment along the 
symmetry axis is related to the principal moments of inertia in a special way. 

1 Introduction 

Let us consider a rigid body B with mass m and centre of mass Oi moving in the 
gravitational field of a point O with mass M, see Figlll We assume that the orbit 
is circular and that it lies in the (x, i/)-plane in the inertial reference frame defined 
by the orthonormal versors {ei, 62, £3} with the origin at O. The principal axes ref- 
erence frame of the body with the origin at Oi is given by the orthonormal versors 
{ai, a2, fls}. We describe the rotational motion of the body with respect to the orbital 
reference frame {s, t, n} with the origin at Oi. Its axes lie along the radius vector of 
the centre of mass of the body, the tangent to the orbit in the orbital plane, and the 
normal to the orbital plane, respectively. 
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Figure 1: A rigid satellite in an orbit around a gravitational centre 

We accept the following convention, see fT]. For a vector q we denote by Q = 
[Qi/ Q2, QsV the associate coordinates in the body frame, i.e., Q, = a, • q, for i = 
1, 2, 3. For two vectors q and p we denote their scalar and vector products hy q ■ p 
and q X p, expressed in terms of their coordinates in the body frame by (Q, P), and 
[Q, P], respectively. Thus we have 

(Q,P) := £Q,.p,. = Q^P = ^.p, 
i=l 
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and 



The equations of the rotational motion of the body can be written in the following 
form 





'Q2P3 


- Q3P2' 




{q X p) 


■ «l" 


Q,P] := 


Q3P1 


-Q1P3 




iqxp) 






Q1P2 


-Q2P1 




_{q X p) 


■ «3. 



dt 



[M,a]+p, 



^N=[N,n], 



(1) 



[S,n-a;oN], 



where M, O. := I~^M, I := diag(A, B, C) are the angular momentum, the angu- 
lar velocity and the inertia tensor of the body, respectively; cuq denotes the orbital 
angular velocity of the centre of mass of the body and P is the torque acting on the 
body. The explicit form of P depends on a particular model. The gravity-gradient 
torque is usually approximated by the following formula 



PG:=3a;|[S,IS], 



where 



w 



_ GM 



and r is the radius of the orbit, see Let us note that in the case of a circular 

Keplerian orbit cuq = w^- Examples of models with cuq 7^ cvk can be found in 



Kepierii 

Eilzit]. 



In this paper we consider the case when, in addition to the gravitational torque, 
also the magnetic torque plays a significant role. Namely, we assume that the grav- 
ity centre (the Earth) is the source of a magnetic field which can be well approxi- 
mated by a magnetic dipole whose axis coincides with e^. Modelling of the mag- 
netic torque Pm is generally difficult because it depends not only on the presence 
of constant magnets located in the satellite, but also on magnetic and conductive 
properties of the material used for its construction, as well as on the presence of 
electronic equipment, for details see IH]. In this paper we assume that the magnetic 
moment of the satellite is induced by the magnetic field of the central body, and, 
moreover, that the body is magnetically symmetric along an axis I fixed in the body. 
Then we have 

Pm:=^(L,N)[L,N], 

where ^ is a parameter depending on the strength of the central magnetic dipole 
and magnetic properties of the body. 
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Thus, we consider the following system 



— M = [M, a] + 3cvi[S, IS] + ^ (L, N) [L, N], 

^N=[N,n], 



dt 



s = [s,n- cuqN]. 



(2) 



It possesses the Jacobi type first integral 

H = i(M,I-iM) - a;o(M,N) + |a;|(S,IS) - ^^(L,N) 

and three geometric first integrals 

H2=(S,S), H3=(N,N), H4=(N,S). 
The above equations can be rewritten in the Hamiltonian form 



(3) 



(4) 



^M, = {M„ H}, ^N, = {N,, H}, ^S, = {S^, H}, 
where the Poisson bracket {•, •} is defined by 



z = 1,2,3, (5) 



3 3 

{M^Mj} = - £ £///cMfc, {M„Nj} = - £ e.j^N^, 
fc=l k=l 

{M„Sj} = - £ e.jkSk, {N,,Nj} = {S„Sj} = {N,,Sj} = 0, 

k=l 



(6) 



where £,y;t is the Levi-Civita symbol. This Poisson bracket is degenerated and the 
three geometric integrals (|4|l are its Casimirs. Their common levels are symplectic 
manifolds [23]. From the geometric interpretation of the vectors N and S it follows 
that, for further study, we can select the following six dimensional symplectic leaf 



= {(M, N, S) G M*^ I (S, S) = 1, (N, N) = 1, (N, S) = 0}, 
which is diffeomorphic to x SO (3, M). 



(7) 



Remark 1. The configuration space of a rigid body whose centre of mass moves in 
a prescribed orbit is SO (3, M) — all possible orientations of the body with respect 
to the orbital frame. Thus the classical phase space of the system is T*SO(3, M) ~ 
X SO(3,M). 
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Remark 2. We can look at system ^ as a Hamiltonian system defined on a nine di- 
mensional Poisson manifold which is s* — the dual to nine dimensional Lie algebra 
5 = (M^ X M^) X so(3, M) (here x denotes the semi-direct product of Lie algebras). 
Then the Poisson bracket defined by ^ is the standard Berezin-Kostant-Kirillov- 
Souriou bracket, and Ai^ is a co-adjoint orbit, see [23]. Here we refer the reader to 
paper [3 ] where the case of a rigid satellite without the influence of magnetic torques 
is considered. 

System (O depends on the parameters p := (A, B, C, cvq, wy^, Li, L2, L3, £,). They 
belong to a set 

7? := {p G X I (L,L) = 1, A<B + C, B < C + A, C<A + B,}, 

whose interior is an eight-dimensional subset of x M"^ (R+ denotes the positive 
real axis). 

It is natural to ask for which p ^ V system ^ or its restriction to admits 
one or two additional first integrals. The high dimensionality of the system and a 
big number of parameters make this problem very difficult. Let us enumerate some 
known facts. 

1. For £, = (the magnetic torque vanishes) the only known completely inte- 
grable case is a spherically symmetric case A = B = C. This case is trivial 
because for a spherically symmetric body the gravitational torque vanishes. 
There is no proof that system (|2)| is non-integrable when ^ = and the body is 
not spherically symmetric. 

2. For ^ = and an axially symmetric body, e.g. A = B, system (O admits one 
additional first integral, namely H5 = M3. There is no proof that this is the 
only situation when system (O possesses one additional first integral. 

3. For A = B = C only the magnetic torque acts on the body. System (O is 
completely integrable and the additional first integrals are H5 = (M, N) and 

= (M, L). In this case the first two equations form a closed subsystem 
which coincides with a special case of the Kirchhoff equations for a rigid body 
in ideal fluid in the integrable case of Clebsh, see Ll^. 

Some limiting cases of system Q when cuq = 0, or cf^ = are worth mentioning 
because they are related to very well known systems. 

Let us consider the case cuq = 0. Now, system ^ describes the rotational motion 
of a rigid body with the mass centre fixed in the external gravity and magnetic fields. 
For ^ = the first and the third equation in (|2)( form a closed subsystem which coin- 
cides with the equations of motion of the completely integrable Brun problem iflll] . 
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see also When cuk = 0, a subsystem of ^ consisting of the first two equations, 
is again a special case of the Kirchhoff equations, see [18]. 

The aim of this paper is to study the integrability of system (O when the body 
is axially symmetric. For this purpose we api^ly the Morales-Ramis theory |25l: l26(1 
which is an extension of the Ziglin theory [3t; 's^. Both theories are based on a 
study of variational equations around a particular non-equilibrium solution of the 
complexified system. We can associate with the variational equations the mon- 
odromy and the differential Galois groups. When the system is integrable, then 
these groups are of a special form and this fact gives a necessary condition for inte- 
grability. To make the paper self-contained, we present basic theoretical facts con- 
cerning the Ziglin and Morales-Ramis theory in the next section. More technical 
material needed in our investigation is presented in the Appendix. We present both 
theories trying to avoid formal language, and we give several examples, which, as 
we hope, helps to understand basic notions of both theories and to popularise them 
in the celestial mechanics community. It is worth mentioning that one of the most 
difficult problems of celestial mechanics — the question about the non-integrability 
of the three-body problem — has been recently solved with the help of these theories, 
see We remark here that H. Poincare lE^ investigated the question 

of integrability of the three-body problem however he assumed that the first inte- 
grals are holomorphic functions of the perturbation parameter (mass of one body). 
Thus, his non-integrability theorems do not assert anything for fixed value of this 
parameter. 

In Section |3l we derive the variational equations along a family of particular so- 
lutions. Our first non-integrability theorem is formulated and proved in Sectional 
We show in this section that the complexified system considered does not possess 
an additional complex meromorphic first integral which is functionally independent 
from the Hamiltonian. The question whether the system does not possess an addi- 
tional real meromorphic first integral is much more difficult. We investigate it in the 
last Section combining the differential Galois approach with the Ziglin argumenta- 
tion 1221] • 



2 Theory 

In this section we describe informally basic facts concerning the Ziglin and Morales- 
Ramis theories. For detailed exposition we refer the reader to lA 25: .21. 
Let us consider a complex dynamical system 

= v{x), teC, X G M", (8) 
where M" is a complex n-dimensional analytic manifold (we can think that M" is 
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just C"). If <p{t) is a non-equilibrium solution of then the maximal analytic con- 
tinuation of (p{t) defines a Riemann surface V with i as a local coordinate. Here it is 
important to distinguish between the abstract Riemann surface V and its image z(r) 
in M". It is crucial when the global geometric language is used. The importance of 
this distinction is discussed in [|24^] . 

Example 1. If (p{t) is given by rational functions oft then V is the Riemann sphere 
CP^ with some points removed (poles of (p{t) ). 

Example 2. If(p{t) is given by elliptic functions with fundamental periods Ti and 
then V is a torus T with some points removed (poles of (p{t)). Moreover, T = C/L, 
where L = {z e C\z = iTi + ;T2, {i, j) G Z^}. 

Together with system ^ we also consider the variational equations 

£, G C". (9) 



= A{t% A{t) = ^(<P(0). 



Let us note that one solution of the above system is known. In fact, if we put 77 
v{(p{t)), then 



dv 



cp{t))v{cp{t)) = A(Ot]. 



(10) 



dt' 3x ^'^M^ ^ dx 
Example 3. Let us assume that system © admits the following invariant set 

n= {{xi,...,Xn) eC"\xi = ■■■ = Xn-l =0}, 

i.e., the right hand sides v{x) = {vi{x),. . . , u„ (x) ) of © are such that Vi (x) = for 
i = 1, ... n — 1 when = for i = 1, . . .n — 1. Then a particular solution (p{t) lies 
on the n-th coordinate axis. Obviously, we have 



dvi 
dXn 



{(p{t))=0, i = l,. 



n 



1. 



Thus, the matrix A{t) has the following block form 

Ait) 



bit) 





ait) 



(11) 



where 
Bit) 



dVi 
dXj 



bit) 



dVn 
dXj 



ait] 



dvn 

dXn 



i,j = l,...,n 



Thus, the first n — 1 variational equations form a closed sub-system of equations 
which are called the Normal Variational Equations (NVEs). 
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The above example shows that the order of © can be reduced by one, at least lo- 
cally. However, these local reductions can be performed consistently over the whole 
r, so we can talk about the NVEs associated with V. For a global definition of the 
NVEs see (37';'?]. Here, just for simplicity, we assume that the coordinates x are 
chosen as in Example |3l Thus, the NVEs have the form 

= B{t)l I = (^1, . . . , 4-1) e C«-i, (12) 
where B{t) is the (n — 1) x (n — 1) upper diagonal sub-matrix of matrix A{t), see 

Remark 3. If system © is Hamiltonian then n is even (n = 2m ) and we have one 
first integral, namely the Hamiltonian of the system. Then we can reduce the order 
of the variational equations by two. Let for our particular solution the value of the 
Hamiltonian be E. Then we can restrict ^ to the level H(x) = E, and we obtain a 
system of 2m — 1 autonomous equations with the same particular solution. Then we 
perform the above-mentioned reduction of the corresponding variational equations 
(of order 2m — 1), and we obtain the NVEs of order 2{m — 1) which are Hamiltonian 
ones. The last statement follows from the Whittaker theorem about isoenergetic 
reduction of order of a Hamiltonian system. 

Remark 4. A typical situation with a Hamiltonian system is the following. For 
the investigated system with Hamiltonian function H{x), x = {q\, pi, . . . , q-m, Pm) £ 
C^"^ there exists an invariant canonical plane T\, e.g., 

n = {{qi,Vl'---'^m,Vm) e C^"' \ qi = Vl = ■ ■ ■ = ^m-l = Pm-l = 0}. 
This implies that 

dH, , an, , ^ „ . , 

— ix) = - — ix) = 0, X G IT, z = l,...m — 1. 
dqi dp, 

Thus, the Hessian of H calculated for x G TT has the following block form 

^""^ - [ hm{x)_ 

where h{x) is a symmetric 2(m — l)x2(m — 1) matrix, and hm{x) is a symmetric 
2x2 matrix. For a particular solution (p{t) G TT the variational equations have the 
form 

E, = ]mH"{(p{t))E„ E G C^'", 

where }m is the symplectic unit (of dimension 2m x 2m ), and the normal variational 
equations are the following 
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Example 4. Let us consider the Hamiltonian system given by the following Hamil- 
tonian function 

H=l{pl + pI) + liql + 4) + (13) 

where e is a parameter and {qi, pi, q2, pi) ^ C"^. The Hamilton's equations for this 
system admit the following particular solution (p{t) = {0,0,q2{t)>P2{t)), where 
q2it) = cn(f,?c), p2(0 = ~ sn{t,k) dn{t,k), k = -s/2/2, and sn, cn, dn denote the 
Jacobi elliptic functions. As this particular solution lies in the {q2, P2) pl^ne, the 
NVEs correspond to variations inqi and pi, so they have the following form 

4 = T], 77 = -eq2{tf£,. (14) 

Note that the above system is a Hamiltonian one. It is generated by the time depen- 
dent Hamiltonian function h = (77^ + eq2{t)'^£,^)/2. 

In the Ziglin and Morales-Ramis theories the concepts of the monodromy group 
and the differential Galois group play fundamental role. In the successive subsec- 
tions we introduce these concepts and give formulations of basic lemmas and theo- 
rems which we used in this paper. 

2.1 Monodromy group 

Let E(i) be the matrix of fundamental solutions of (|9|l defined in a neighbourhood 
of tQ G C, i.e., columns of are n linear independent solutions of and let y 
be a closed path (with the base point at to) on the complex time plane. An analytic 
continuation of E(^) along y gives rise to a new matrix of fundamental solutions Z{t) 
in a neighbourhood of to which does not necessarily coincide with 'El{t). However, 
the solutions of a linear system form an n dimensional linear space, so we have 
E(f) = E(f)My, for a certain nonsingular matrix My G GL(n, C) which is called the 
monodromy matrix. 

Example 5. The system 



d 


'^1 


1 


'0 






di 


£,2 


~ i2 


-1 t 




£,2 



has two linearly independent solutions 

^(i) = (i,l)r, and ^^^^ = {tlnt,l+'^nty. 

After continuation along a loop y encircling t = once, the solution ^(^^ is un- 
changed. However, the second solution changes into 
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and thus we have 



tint 
1 + lnt 



E(f)M. 



r 



't tint ' 




'1 2m 


1 1 + lnt 




1 



Example 6. Let us consider the following system 

where C is a constant matrix. Let ybe a loop encircling once t 
Then the monodromy matrix is given by 



counterclockwise. 



My = exp[27riC]. 

The monodromy matrix My does not depend on a particular choice of y. If the 
path cr can be obtained by a continuous deformation of the path y, then M^ = My. 
We denote by [y] the set of all paths which can be obtained by continuous defor- 
mations of y, and it is called the homotopy class of path y. Thus, the monodromy 
matrix My depends on the homotopy class of path y. If we have two paths cr and y 
by their product r = a ■ y we understand the path t obtained in the following way: 
first we go along y, then along cr. One can show that this defines properly a product 
of homotopy classes, i.e., [t] = [cr] • [y] := [cr • y]. We can also define the inverse 
y~^ of the path y: we go along y in the opposite direction. Again we have a correct 
definition [y]~^ := [y~^]- In this way the homotopy classes form a group which is 
called the first homotopy group of a Riemann surface (walking on the complex time 
plane t, in fact we make loops on V because t parametrises the surface P). We denote 
it by 7ri(r, ^o)- 

Remark 5. If we change the base point tQ of the paths, then, instead of the matrices 
My, we obtain CMyC^^, where C is a certain nonsingular matrix (the same for all 
paths). It means that the homotopy groups at all points to are isomorphic. 



All the monodromy matrices form a group Ai with respect to matrix multipli- 
cation which is a subgroup of GL(n, C). From the definition of monodromy we 
have Mff.y = MyMa, so Mj^,.^] = M[^]M[o.]. In the same way M[^]-i = M^|. In 

other words, the monodromy matrices form an anti-representation of 7X\{V,tQ) in 
GL(n,C). 



Remark 6. If system Q is Hamiltonian, then the variational system © is also a 
Hamiltonian one, and the monodromy group is a subgroup of the symplectic group 
Sp(2m, C) where 2m = n. If we consider the NVEs for a Hamiltonian system as 
it was described in Remark 13 then the monodromy group of these equations is 
contained in Sp(2(m — 1), C). 
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2.2 Basic lemma of the Ziglin theory 

Let us assume that F{x) is a holomorphic first integral of ©. The Taylor expansion 
of F{(p{t) + E,) has the form 

F((p(i) + ^) = F((p(i)) +F^(i,^) + • • • , (15) 

where Tm{t, £,) is a homogeneous polynomial (with respect to the coordinates of E) 
of degree m > 0. It is easy to show that Fm{t, K) is a first integral of the variational 
equations We called Fni{t,E,) the leading term of the first integral. When the 
first integral F{x) is a meromorphic function, then it can be represented as a ratio 
P{x)/Q{x) of two holomorphic functions P{x) and Q{x). If Pm{t,E,) is the leading 
term of P{x) and Q]^{t, E) is the leading term of Q{x), then by the leading term of 
f (x) we understand Pm{t,£,)/ Qk{t, E,), and it is a first integral of equations ^ which 
is rational with respect to ^. 

An analytic continuation of solutions of ^ along a closed path y transforms 
initial conditions for these solutions to other points in the following way At Iq we 
start from ^q- For small t we move along y and goes to E,{t) = Z{t)£,Q. After 
continuation, we return to a neighbourhood of to, but now our point is moved to 
Z{t)£,Q, and thus at the end of the path at we obtain the point 

as E(io) is the identity Thus we have the following map 

(^o,^o) ^ (fo,My^o)- 

r 

It is important to notice here that t^, as well as E,q, are arbitrary. 

Let Frn{t, E) he a first integral of © and let = Fm{to,Eo). A first integral does 
not change its value when we make an analytic continuation. Thus taking the loop 
y we have 

As tQ, E,Q and y are arbitrary we have 

Fm{t,E,)=Fm{t,Mr^), (16) 

for all My G A^. In other words, Fm{t,E) is invariant with respect to the natural 
action of the monodromy group. A non-constant function satisfying the above con- 
dition is called a first integral (or an invariant) of the monodromy group (polynomial 
(rational) if Fm is a polynomial (rational) function of the coordinates of £,). We can 
repeat all the above considerations for the normal variational equations. The condi- 
tion ((T6)l is restrictive. When the monodromy group of the NVEs is 'big', then it can 
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happen that there is no non-constant polynomial (rational) invariant, and this fact 
implies that system ^ does not have a holomorphic (meromorphic) first integral. 

The following lemma formulated by Ziglin gives the necessary condition for in- 
tegrability, see Proposition on p. 183 in |37H and Proposition on p. 4 in jsitl . 

Lemma 1. If system ^ possesses a meromorphic first integral defined in a neighbourhood 
U C M", such that the fundamental group of V is generated by loops lying in U, then the 
monodromy group Ai of the normal variational equations has a rational first integral. 

Remark 7. The reason why in the above Lemma the necessary condition for inte- 
grability cannot be formulated (or, rather, it is more difficult to formulate) in terms 
of the monodromy group of the full variational equations is the following. The 
monodromy group of ^ always possesses one polynomial invariant. Let us ex- 
plain why. As it was mentioned, for equations ^ we know one particular solution 
r\ = v{(p{t)), see (ITOt . IfZ{t) is the fundamental matrix of then we can find a 
vector c G C" such that r\ = r]{t) = Z{t)c. Let us assume for simplicity that the 
solution (p(t) is single-valued. Thus the continuation ofr\{t) along an arbitrary path 
7 does not change it, and we have thatr\{t) = 'E.{t)c = Z{t)MyC = Z{t)c. It follows 
that MyC = c, i.e., the vector c is an eigenvector of all monodromy matrices and it 
corresponds to an eigenvalue 1. Thus, in appropriate coordinates x = (xi, . . . , x„), 
the monodromy matrices M can be put simultaneously into the following form 



M 



1 

m M 



where m, M are (n — 1) x 1 and {n — 1) x (n — 1) matrices, respectively. But now 
the linear polynomial f{x) = Xi is an invariant of the monodromy group. 



2.3 Differential Galois group 

Let us assume that the entries of the matrix A{t) of the linear system ^ are rational 
functions of t. We know that solutions of linear equations with rational coefficients 
are not necessarily rational, however, we can ask whether a given linear equation or 
a system of linear equations is solvable in terms of 'known' functions. This question 
was investigated at the end of the nineteenth and at the beginning of the twenti- 
eth century by Picard, Vessiot and others. Later on, thank to works of Kolchin, the 
Picard-Vessiot theory was considerably developed to what is now called the differ- 
ential Galois theory. For a general introduction to this theory see i2j;E3iSl23l- 

Through this subsection our leading example is a linear second order differential 
equation with rational coefficients 

w" + pw' + qw = 0, p,qeC{t), ' = (17) 
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In what follows we keep algebraic notation, e.g., by C[t] we denote the ring of poly- 
nomials of one variable t, C(i) is the field of rational functions, etc. Here we consider 
the field C{t) as a differential field, i.e., a field with distinguished differentiation. 
Note that in our case all elements a G C{t) such that a' = are just constant, i.e., we 
have fl' = G C. Thus such elements form a field — the field of constants. 

Remark 8. In the most general case we meet in applications, the coefficients of © 
are meromorphic functions defined on a Riemann surface V, which is usually de- 
noted by Ai{r). Meromorphic functions on V form a field. It is a differential field if 
equipped with ordinary differentiation. 

The field C(^) can be extended to a larger differential field K such that it will 
contain all solutions of equation (fTTjl . The smallest differential field K containing n 
linearly independent solutions of © is called the Picard-Vessiot extension of C{t) 
(additionally we need the field of constants of K to be C). 

Remark 9. The Picard-Vessiot extension for equation ((T7)l can be constructed in the 
following way. We take two linearly independent solutions £, and r\ of ((T7)l (we 
know that such solutions exist). Then, as K we take all rational functions of five 
variables {t, ri, ri'), i.e., K = C{t, ri, ri'). 

Remark 10. In the case considered (a system of complex linear equations with ra- 
tional coefficients) the existence of the Picard-Vessiot extension follows from the 
Cauchy existence theorem. In abstract settings, i.e. when we consider a differen- 
tial equation with coefficients in an abstract differential field, the existence of the 
Picard-Vessiot extension is a non-trivial fact, see e.g. [22]. 

Now, it is necessary to define what we understand by 'known' functions. Infor- 
mally, these are rational and algebraic functions, their integrals and exponential of 
their integrals. More precisely, we say that a solution 77 of ((TTjl is: 

1. algebraic over C(^) if t] satisfies a polynomial equation with coefficients in C(f), 

2. primitive over C(^) if rj' G C(^), i.e., if 77 = fa, for certain a G C{t), 

3. exponential over C{t) if vj' /rj G C{t), i.e., if 77 = exp / a, for certain a G C{t). 

We say that a differential field L is a Liouvillian extension of C(i) if it can be obtained 
by successive extensions 

C{t) = KoCKiC-- - CKm = L, 

such that Ki = Ki_i{rii) with 77/ either algebraic, primitive or exponential over Ki_i. 
Our vague notion 'known' functions means Liouvillian functions. We say that Q is 
solvable if for it the Picard-Vessiot extension is a Liouvillian extension. 
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Remark 11. All elementary functions, like e , log t, trigonometric functions, are 
Liouvillian, but special functions like Bessel or Airy functions are not Liouvillian. 

Example 7. The equation 

4tw" + 2w' -w = 0, 

has two linearly independent solutions Wi = exp[y/t] andw2 = exp[—\/t]. Both of 
them are Liouvillian. 

How can we check if solutions of a given equation are Liouvillian? For this 
purpose we need to check properties of the differential Galois group of the equation. 
This group can be defined as follows. For the Picard-Vessiot extension K D C{t) we 
consider all automorphisms of K (i.e. invertible transformations of K preserving 
field operations) which commute with differentiation. An automorphism g : K ^ K 
commutes with differentiation if g{a') = {g{a))' , for all a ^ K. We denote by A the 
set of all such automorphisms. Let us note that automorphisms A form a group. 
The differential Galois group Q of extension K D C(f), is, by definition, a subgroup 
of A such that it contains all automorphisms g which do not change elements of 
C{t), i.e., for g G ^ we have g{a) = a for all a G <C{t). 

Remark 12. It seems that the definition of the differential Galois group is abstract 
and that it is difficult to work with it. However, from this definition we can de- 
duce that it can be considered as a subgroup of invertible matrices. Let Q be the 
differential Galois group of equation ((T7)l and letg^Q. Then we have 

= ^(0) = g{^" + v^' + = g{^") + g{v)g{^') + g{^)gip)' 

but g commutes with differentiation so g{w") = {g{w))" , g{'w') = {g{w))',and, 
moreover, g{p) = p, g{q) = q because p,q G C{t). Thus we have 

(gH)" + pigH)' + igM = 0. 

In other words, ifw is a solution of equation ((T7|l then g{w) is also its solution. Thus, 
if B, and r\ are linearly independent solutions of (fTTjl , then 

g{^) = gn^ + ^21^7/ g{^) = gi2^ + gim, 

and 



g 



( 










'gw 


gl2 








I' vl 




g21 


g22_ 



Hence, we can associate with an element g of the differential Galois group Q an 
invertible matrix [gij], and thus we can consider Q a subgroup of GL (2, C) . If instead 
of the solutions E, and r\ we take other two linearly independent solutions, then all 
matrices [gij] are changed by the same similarity transformation. 
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The construction presented in the above remark can be easily generalised to a 
linear differential equation of an arbitrary order and to a system of linear equations. 
Thus we can treat the differential Galois group as a subgroup of GL(n, C). Let us 
list basic facts about the differential Galois group 

1. If g{a) = a for all g ^ G, then a G C{t). 

2. Group Q is an algebraic subgroup of GL(n, C). Thus it has a unique connected 
component which contains the identity, and which is a normal subgroup of 
finite index. 

3. Every solution of the differential equation is Liouvillian if and only if Q'^ con- 
jugates to a subgroup of the triangular group. This is the Lie-Kolchin theorem. 

For proofs and details we refer the reader to the cited references. 
2.4 Basic theorem of the Morales-Ramis theory 

For a given linear system of linear differential equations we can determine the mon- 
odromy group Ai and the differential Galois group Q. From the description given 
above it follows that both these groups are related. In fact, we have C ^. In 
other words, the differential Galois group Q is 'bigger' then the monodromy group 
M. 

Example 8. For the Airy equation x = tx the monodromy group is trivial, i.e., it 
contains only one element — the identity matrix, while its differential Galois group 
is SL (2, C) . For a proof see e.g. [1 6 ]. 

Remark 13. It should be mentioned that the determination of the monodromy 
group is a difficult task, and this groups is known only for a very limited num- 
ber of equations. What concerns the determination of the differential Galois group 
we are in much better situation. There exist algorithms (the Kovacic algorithm [17]) 
which allow to determine this group for an arbitrary second order linear differential 
equation with rational coefficients (see the Appendix for additional references). 

The fact that Ai C G suggests the use of Q instead of Ai to formulate a necessary 
condition for non-tntegrability. If system ^ possesses a meromorphic first integral, 
then ^ has a first integral and this fact imposes a restriction on its differential Galois 
group Q, as it imposes restrictions on its monodromy group Ai. In fact, we have a 
lemma which is analogous to Lemma [H 

Lemma 2. If system ^ possesses a meromorphic first integral defined in a neighbourhood 
UcMofV, then the differential Galois group Q of the NVEs has a rational first integral. 
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The above lemma is a variant of Lemma III.1.13 from [2], see also Lemma 4.6 in 
lE^ . For proof and details see Chapter III of 10] • 

The differential Galois theory gives a powerful tool to the study of integrabil- 
ity of Hamiltonian systems. The Morales-Ramis theory is formulated in the most 
exhaustive form in book fzs'] and papers [26]. It gives a necessary condition of in- 
tegrability of a Hamiltonian system for which we know a non-equilibrium solution. 
The main theorem is the following 

Theorem 1. Assume that the Hamiltonian system is integrable in the Liouville sense in 
a neighbourhood of a particular solution. Then the identity component of the differential 
Galois group of the NVEs is Abelian. 

3 Particular solutions and variational equations 

From now on we consider (|2|l as a complex system, i.e., we assume that (M, N, S) G 
and t G C. Without loss of generality, choosing appropriately the unit of time 
and length, we can put cuk = <^o = 1 ^rid A = 1. According to our knowledge, 
for an arbitrary L, equations ^ do not admit a particular solution. However, if we 
assume that L coincides with one of the principal axes, e.g., L = [0, 0, 1] ^, then one 
can find particular solutions. In fact, in this case the following manifold 

M={{M, N, S) G C'^ I M2 = Ms = N2 = N3 = Si = 0, Ni = 1}, (18) 

is invariant with respect to the flow generated by system (|2|l. Solutions lying on J\f 
describe the planar rotations of the satellite when its third axis is permanently in 
the orbital plane and its first axis is perpendicular to the orbital plane. Moreover, 
we can easily find an analytic form of the solutions of (|2)| describing this motion. In 
fact, system ^ restricted to J\f has the form 

Ml =3(C-B)S2S3, S2 = (ai-l)S3, S3 = -(ai-l)S2, (19) 
and it possesses two first integrals 

H|A^ = ^M? - Ml + ^ (bsI + CSl) , H2|^ = + si (20) 

We can introduce on the level H2 = 1 local coordinate (p such that 

S2 = — cos (j) and S3 = sin 4>. 

Then system ((T9)l reads 

Ml = -3(C-B)sin(/jcos(^, = Mi - 1. (21) 
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Thus, we have 

(jb = — 3(C — B) sin (p, (p = 24>. (22) 

Solving the above equation we obtain an one parameter family 0(i, fc) of the solu- 
tions of ^ expressed in terms of the Jacobi elliptic functions. Let us define 



cv=^3\C-B\. (23) 
Then the explicit form of the solutions is given by 

Mi{t,k) = 1 + cvkcn{cvt,k), (24) 



and for C > B 



for C < B we have 



Siitfk) = — dn{cvt,k), 
S3{t,k) = ksn{ajt,k); 

Siit'k) = ksn{cvt,k), 
S3{t,k) = dn{cvt,k), 



(25) 



(26) 



where 



e (0, 1), (27) 
and E is the value of the energy integral for equation (|22)l. i.e., 

r ^ ■ 2 2 

E = -(p — cu cos (p. 

Let us note that for the above solutions we have 

H{0{t,k)) = \w\^ + 1(3B - 1) := h{k). (28) 

From the above formulae it follows that the particular solutions given above are 
single-valued, meromorphic, and double periodic with periods 

T{k) = -K{k), T{k) = -iK'{k), 
cv cv 

where K{k) is the complete elliptic integral of the first kind with modulus k, K'{k) := 
K{k'), and k' := \/l — k^. In each period cell they have four simple poles at: 

Ti(fc) = Ir(fc) + iT'(fc), T2(fc) = Ti(fc) + Ir(fc), 

^ T T (29) 

T3(fc) = T2(fc) + -T{k), r,{k) = rs{k) - -T'{k) mod (T(fc), T'{k)). 
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Thus, the Riemann surfaces associated with the particular solutions ^{t,k) are 
tori with four points: S/(fc) = 0(T;(fc),fc), / = 1,2, 3,4 removed. In with coordi- 
nates (Ml, S2, S3) these Riemann surfaces are intersections of two quadrics 



For < k < 1 the four points S/ {k) correspond to four points of intersections of the 
above quadrics at infinity. 

As our aim is to investigate the case when the satellite is symmetric, we assume 
that A = B = 1. For a symmetric satellite, we have one more first integral, namely 
H5 = M3. This first integral is connected with the existence of an one parameter 
symmetry of the system. Equations ^ are invariant (for the prescribed choice of 
L and the symmetry axis) with respect to an action of group SO (2, M). Simply, the 
principal axes perpendicular to the symmetry axis of the body can be chosen arbi- 
trarily. Thanks to that, we can reduce the number of degrees of freedom by one. 
Thus, the reduced system is Hamiltonian with two degrees of freedom and it de- 
pends parametrically on the value of the chosen level of H5. 

Further calculations can be performed in the (M, N, S) coordinates in the same 
way as it was done in ||^. Here we perform them in canonical coordinates on 
This approach allows to deduce the normal variational equations in an elementary 
way. Appropriate canonical variables on can be chosen in the following way. 
We parametrise the orientation of the principal axes of the body with respect to 
the orbital reference frame by the Euler angles {q\,q2r ^3) of the type 3-1-3, and we 
take them as generalised coordinates. Then generalised momenta conjugated to 
('?i/'?2/^?3) are given by 





(30) 



p = KM, 



K 



sm sm q2 cos sm q2 cos q2 
cos q3 — sin qs 



(31) 



1 



Moreover, we have 



sin qs sin q2 — sin qj cos q2 sin qi + cos qs cos qi 

N= cosi^3sini^2 S= — cos i^3 cos siri — sin i^3 cos t^i 

cosi^2 sini^2sini?i 



(32) 



In the introduced canonical coordinates the Hamiltonian ^ reads 




(33) 
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As we can see, is a cyclic coordinate and p3 = M3 is a first integral. Thus, con- 
sidering p3 as an additional parameter, H defines a Hamiltonian system with two 
degrees of freedom with x = {q\, q2, pi, pi) as canonical coordinates. As our partic- 
ular solutions lie on the level M3 = 0, we investigate this system for p3 = , i.e, we 
consider the Hamiltonian system given by the following Hamiltonian 

H = + hl-pi + ^{C- 1) sin2 q^ sin^ q2 - h cos^ q2. (34) 

Now, the invariant manifold J\f corresponds to the canonical plane q2 = 7i/2, p2 = 
0, on which canonical equations generated by H have the form 

q\ = p\ — l, p\ = —3{C — l)su\qicosq\. (35) 

Comparing them with equations (|2T|l we see that pi = Mi and qi = 4> (note that 
this fact follows from the definition of A/", formulae (|3T|l . ((32)l and the fact that on J\f 
we have (^3 = 7i/2). Thus, the explicit form of the particular solutions x = x{t,k) = 
(qi (t, k), n/2, pi {t, k), 0) is given by 



and for C > 1 



and for C < 1 



pi{t,k) = 1 + cukcn{cvt,k), (36) 

cos qi {t, k) = dn(cuf, k), 
sm.qi{t,k) = ksn{wt,k), 



cosqi{t,k) = -ksn{cut,k), ^^^^ 
smqi{t,k) = dn(cvt,k). 



We note here that for a symmetric satellite we have 



3|C- 1| G (0,\/3). 



Cf = 

The variational equations along the particular solution x{t,k) have the following 
form 

Qi = Pi, Pi =3(l-C)cos(2^?i(i,fc))Qi, (39) 
Q2 = P2, P2=[^- pi{t,kf + 3(C - 1) sin2^^i(^,fc)]Q2. (40) 

As the particular solutions lie in the plane {q2 = tc/2, p2 = 0}, the NVEs correspond 
to the subsystem (l40t which can be written as a second order linear equation of the 
form 

Q + fl(i,fc)Q = 0, Q = Q2, (41) 
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where 

a{t,k) 



[l + kcucn{cut,k)f + cu^drP{cut,k) - E, for C < 1, 
[l + kcvcn(cvt,k)f - cv'^k'^siP{cvt,k) - £, for C> 1. 



(42) 



Remark 14. Let us notice that for equation (|4T|I the differential Galois group is a 
subgroup of SL(2, C) . It is always the case when a second order linear differential 
equation does not contain a term proportional to the first derivative. 



Remark 15. Here we underline that the obtained NVE (|4T|I is the reduced normal 
variational equation for (O when A = B = 1 and L = (0,0,1) derived for so- 
lution (|24|) - ((26)l . We just performed a symplectic reduction as in [3], but for this 
purpose we use local canonical coordinates. 



Equation (|4T|I is defined on Vj^. In order to use the differential Galois theory ef- 
ficiently, it is crucial to transform the investigated equation into an equation with 
rational coefficients. In our case we can do this making the following transforma- 
tion 

t — > z := kcn{cvt,k). 

Then the NVE has the form 



(43) 



where 



p(z) 



q{z) 



Q" + p{z)Q' + q{z)Q = 0, 
z{-l + 2{k^-z^)) 
' + (1 + cvzf + aP- [z^ + k'^' 

aP (fc2 - z2) (^z2 + fc'2 j 

+ (1 + cuzf - aP (fc2 - z2) 



d_ 

Tz 



(44) 



CXP (fc2 



z2 + k! 



for C < 1, 



for C> 1. 



(45) 



Equation (|44|l is Fuchsian (see Appendix) and it has five regular singular points over 
CP-'^, namely zi^2 = Z3,4 = and Z5 = 00. 

Remark 16. Our transformation ()43)l is a double covering 

cpi — ^ C — > 



The differential Galois groups of equation ()4T|l and equation (@4|l are different, how- 
ever these groups have the same identity components, see ^,25]. 



20 



Changing the dependent variable 



p(s) ds 



Q = Wexp 

we transform (|44|l to the reduced form 
W" = r(z)W, r(z) 
The rational coefficient r(z) has the following simple fraction expansion 



(46) 



(47) 



riz) = t 



+ 



(z - ZiY Z - Zi 



(48) 



with coefficients 



Cl\ = Cl2 — ^3 — 



3_ 
16' 



and for C < 1 



, 3cv^(3 + 4k^)+8{l-^ + 2kcv) , 

bl = TTl 9 / ^2 



h = i 



aP-{l2k'^ + 1) + 8(^ - 1 + lik'w) 
lek'cv^ 



-h + 



cu 



h* 



(49) 
(50) 



where * denotes the complex conjugation. For C > 1 the coefficients hi are the 
following 



h = i 



cv'^jUk^ + 1) + 8(1 - ^ + 2kw) 
16kcv^ ' 
. 3cv^{3 + Ak'^) + 8(^ - 1 + lik'w) 



hi 



cu 



fc4 = hi. 



lek'w^ 

The Laurent expansion of r(z) at infinity in both cases has the same form 



(51) 
(52) 



r(z) 



+ 



1 



(53) 



Remark 17. Transformation (|46|l changes the differential Galois group. For equa- 
tion (|47|) Q is a subgroup of SL (2, C) but for equation (|44)l Q is not a subgroup of 
SL(2, C). Generally, when the coefficients p{z) and q{z) in (|44|) are arbitrary ratio- 
nal functions, transformation (l46b changes also the identity component of Q, e.g. 
of equation (|44)l can be non-Abelian but for the transformed equation (|47)l Q can be 
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Abelian. However, if of equation (|44)l is solvable then of equation (|47)l has the 
same property. In our case transformation (l46t has the following form 

W=[{k^-z^){k'^ + z^)]'/^Q 

and thus it does not change the identity component of the differential Galois group 
of equation (|44|l. This is not accidental. In the time parametrisation the NVE has 
the form (|4T|l and its differential Galois group is contained in SL(2, C). Then we 
make transformation (l43t which is a finite covering, and thus it does not change the 
identity component of the differential Galois group, see Proposition 4.7 in [4 ]. Then, 
by Lemma 4.24 from [4] transformation ()46)l has the form W = RQ, where is a 
rational function for an integer n. 



4 Complex non-integrability 

First, we investigate local monodromy of equation (|47)l at infinity. In many cases it 
simplifies proofs considerably 

Lemma 3. Let us assume that C ^ 1 and 2£, ^ 3(1 — C). Then the local monodromy of 
equation ()47)l at infinity is 

"1 27ril 



Mr 



1 



Proof. We prove the lemma for C < 1. For C > 1 the proof is similar. First we 
change the dependent variable z = 1 /C. This change moves z = oo to C = and 
transforms (|47|) to the form 

W + \W - ^r (l) W = 0. (54) 



Moreover, we have 



1 (1\ 2 
rr 



^4- y^j - ^2+0{r^), (55) 

and thus the indicial equation (see [36, Chapter X]) reads 

p(p- 1) +2p-2 = 0. (56) 

Hence, exponents at C = are p_ = — 2 and p+ = 1. Their difference m = p-|- — 
p_ = 3 is an integer, and thus, in a neighbourhood of C = one solution of (|54|) has 
the form 

00 

Wi(C) = e^f{Q, /(C) = 1 + I/fcC'S (57) 

k= 



11 



where the series defining /(C) is convergent in the considered region Il36l] . The sec- 
ond solution, independent of Wi (C), is defined by the integral 



W2(c) = Wi(0 



'-ds 



Wi(C) 



-m—l 



ds 



f{s) 



2' 



(58) 



Let us denote 



1 + 1 gkC 

k=l 



Then, from ((58)l it follows that the solution W2(C) can be written in the form 

W2{C)=gmWi{C)\nC + e-V{C), (59) 

where V{C) is holomorphic in a neighbourhood of C = 0. The form of local mon- 
odromy depends on whether a logarithmic term is present or not in the solution. 
To check if it is present in our case, we have to calculate if ^ 0. It can be easily 
shown that 

^3 = -2(2/3- 3/1/2 +/3). 

The coefficients /, i=l,2,3 of the expansion ()57)l can be computed directly (see e.g. 
and they are the following 



h = 



tv 



2(4fc2 



l)+2(2-^) 
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a;2(108fc2 -47) +2(9 - 7^) 



360 a;3 



(60) 
(61) 



One can check that 



g3 



2^) 



Thus, if cv^ 7^ 2£, the logarithmic term in the solution W2(C) is present. Note that for 
C < 1 the condition cu^ 7^ 2£, is equivalent to 2E, 7^ 3(1 — C). Now, let us consider a 
small loop y encircling the singular point C = counterclockwise. The continuation 
of the matrix of the fundamental solutions along this loop (under the assumption 
that cu^ 7^ 2E,) gives rise to the triangular monodromy matrix 



Wi(C) W2(C) 



Wi(C) W2(C) 





'1 


2m 







1 



This ends the proof. 



(62) 



□ 



In the next lemma we show that for almost all values of the parameters equa- 
tion (|47)l is not reducible, i.e., for it case 1 in LemmalUdoes not occur. 
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Lemma 4. For C 7^ 1 and G (0, 1) equation ()47)l zs not reducible except for the case when 



3 2 
^ = -(l-C), flnrf '^^ = 2k^-l ' ^^^^ 

Proof. To prove our Lemma we apply directly the first case of the Kovacic algorithm 
(see Appendix). First we consider the case C < 1. All finite poles of r(z) and infinity 
are of the second order. Using the coefficients flf, i = 1, ... ,4 given by (l52t and the 
expansion (|53)| we obtain 

Aoo = 3, (64) 

Eoo = {-l,2}. (65) 

We proceed to the Second Step. From the Cartesian product E = £00 x nf=i we 
select these elements e = {Coo, Si, e2, £3, £4) G £ for which 



d{e) = 1 - 




where No denotes the set of non-negative integers. In our case there exist seven 
elements of £ satisfying this condition 



,(1) = 1 


1111] 

'4'4'4'4J 


d{e^'^) 


= 1, 


= 1 


113 3' 
'4'4'4'4J 


d{e^'^) 


= 0, 


= 1 


3 3 11] 
'4'4'4'4J 


d{e^'^) 


= 0, 


= 1 


^13 3 1] 
'4'4'4'4J 


d{e^'^) 


= 0, 




3 113] 
'4'4'4'4J 


d{e^'^) 


= 0, 


,(6) = 1 


13 13' 
'4'4'4'4J 




= 0, 


,(7) = 1 


3 13 1' 
'4'4'4'4J 


d{e^'^) 


= 0. 



Ai = A2 = A3 = A4 



1 

2' 



and thus 



1 3 
4'4 [' 
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Now we pass to the Third Step of the Kovacic algorithm. For each element e G E 
such that d{e) G Nq, we construct a rational function 

w = w{e) = £ (67) 

i=l 2 

Then we check if there exists a monic polynomial P G C[z] of degree d{e) satisfying 
the equation 

P" + 2wP' + {w' + w^- r)P = 0. (68) 

If we find such polynomial, then equation (|47)l has an exponential solution W = 
Pexp / w. 

For e^'^^ we have d(e^^^) = 1, thus we take P = z + g and, substituting P to 
equation (l68b , we obtain the following algebraic system 



g 



w^{k^ -1) +£,-1 =0, -2ga)-fc^ 0)^ + ^-1 = 0, cv{gw + 1) = 0. (69) 



2 ,,.2 



We note that cu ^ 0, and thus this system has one solution 

1.1 2 2 



For e^''^ i = 2, ... ,7 we have to find a monic polynomial of degree zero satisfying 
((68)l , so we put P = 1. 

For e(2) equation (|68|) yields 

+ 1 - ^ = 0, cu = 0, (70) 

but cf 7^ 0, so there is no solution of the above equations. 
We have the same situation for e^"^) when (l68t gives 

^ - 1 = 0, a; = 0. (71) 

For e*^'") with m = 4, . . . , 7 we obtain two equations of the form 

w^{2k{k T ik') - 3) + 4(^ - 1) = 0, a;(a;(fc ^ ik') ± 2) = 0, (72) 

where the choice of signs depends on m. The second equation cannot be satisfied by 
a real cv ^ Q and A: G (0, 1). This finishes the proof for C < 1. The proof for C > 1 is 
similar. □ 

Combining the above two lemmas we have. 
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Lemma S.IfC^l and 2E, 7^ 3(1 — C), then for k G (0, 1) the differential Galois group 
Gof^is SL(2,C). 

Proof. In fact, under the given assumptions Q cannot be a triangular subgroup of 
SL(2, C) by LemmaHl Under the same assumptions, by Lemma|3l we know that 
contains a non-diagonalisable triangular matrix Moo- Thus case 2 in Lemma|9lcannot 
occur as in this case Q'^ is diagonal. By the same reason case 3 in Lemma |9l cannot 
occur as for a finite group Q the identity component consists of the identity Thus, 
wehave^ = ^0 = SL(2,C). □ 

As SL(2, C) is not Abelian, we have, as a direct consequence of the above lemma, 
the following. 

Lemma 6. If C 7^ 1 and 2E, ^ 3(1 — C), then for k G (0,1) the complexified Hamilto- 
nian system given by (|34)l does not admit an additional complex meromorphic first integral 
functionally independent together with H in a neighbourhood of phase curve F^. 

However, as we mentioned the Hamiltonian system given by (|34)l is a subsystem 
of (|2|l, thus as a corollary we have the following theorem. 

Theorem 2. If C ^ 1, A = B = 1,L = (0,0, 1) and 21 ^ 3(1 - C), then for k G (0, 1) 
the complexified system ^ considered on M.^ does not admit an additional complex mero- 
morphic first integral functionally independent together with H and H5 in a neighbourhood 
of the phase curve r^. 

In the above Theorem the case 2E, = 3(1 — C) is excluded. One can suspect 
that for these values of the parameters our system is integrable. Indeed, the lemma 
below shows that our suspicions are well justified because a necessary condition for 
the integrability is satisfied. 

Lemma 7. If2£, = 3(1 — C) then for all G (0, 1) the identity component of the differential 
Galois group of ()47)l is Abelian. 

Proof. We consider the case C < 1. The proof for the case C > 1 is similar. 

By LemmaHl we know that for 2£, = 3(1 — C) equation (|47|) is reducible only 
when cu^ = 2/ {2k^ ~ !)• As all exponents are rational and equation (|47)l is Fuchsian, 
in this case ^ is a proper subgroup of the triangular group so is Abelian. 

For cv^ ^ 2/(2fc2 - 1) 

we show that case 2 of Lemma IH occurs. To this end we 
apply the Kovacic algorithm for this case. Now sets Ei, E2, E3, E4 and Eoo have the 
following forms 

El = E2 = E3 = E4 = {1,2,3}, Eoo = {-4,2,8}. (73) 
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We have to find at least one monic polynomial P G C[z] of degree 

d(e)=2-if^'oo+£^'«), (74) 



i=l 



satisfying the differential equation 

P'" + 3wP" + {3w^ + 3w' - 4r)P' + {w" + 3ww' + - Arw -2r')P = 0, (75) 
where 

^ = He) = 1 £ (76) 

i=l 

We choose e = (—4, 1, 1, 1, 1). Then d{e) = 2, and we look for a polynomial of the 
second degree 

P{z)=z^+giz + g2, (77) 

satisfying (jTSjl. Substituting (|77|l and (|76)l into (|75)) we obtain the following system 
determining gi and g2 

[(2fc2 - 3)a;2 + 4(^ - l)]gi - 4a;g2 = 0, 3cvgi + 2cv^g2 + 2{k^ + 1 - ^) = 0, 
cv{cvgi + 2) = 0. 

(78) 

If <^ = cu^/2 then the above system has the following solution 

2 (l-2fc2)a;2 + 4 

□ 



5 Real non-integrability 

On J\f system (O has four equilibria 

s± = (1,0,0,0,0,1,0,^1,0), M± = (1,0,0,0,0, 1,0,0, ±1). 

These equilibria correspond to a fixed position of the satellite in the orbital frame. 
For s± the symmetry axis is parallel to the radius vector of the centre of mass of the 
satellite, and for u± the symmetry axis lies in the orbital frame and it is perpendicu- 
lar to the radius vector. 
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7C 2jt 7C 2jt 

Figure 2: Phase portrait of planar oscillations in the {4>,Mi) plane. The left panel 
corresponds to the case when C > 1, and the right panel corresponds to the case 
when C < 1 . 

Let us restrict system ^ to the invariant manifold TT = M H Ai^. Then, for the 
restricted system the equilibrium points u± are hyperbolic if C > 1; if C < 1, then 
s± are hyperbolic. See Figure |2l 

We restrict further discussion to the case C > 1. For real t and < k < 1 the 
solution 0{t,k) defined by (|24)l - ((26)l corresponds to closed phase curves around 
the stable point s+. Closed real phase curves around s_ are given parametrically by 
<I>{t + T' (k) /2), t G R. Let Vi be the phase curve corresponding to the solution given 
by ((24|l - (l26t with k = 1, i.e., 0(t, 1). Then Fi contains four components which are 
real phase curves corresponding to real solutions heteroclinic to u±. Their union is 
Re Fi, and by D. we denote the closure of Re Fi. 

Lemma 8. Let us assume that C ^ 1. Then for an arbitrary complex neighbourhood U C FT 
of D. there exists e > 0, such that for < |fc — 1| < e the fundamental group 7ri(Fjt) of 
phase curve V^- is generated by loops lying in U. 

Proof. The periods T{k) and T'{k) of the solution 0(i, fc) are primitive and at the 
same time they are the minimal real and imaginary periods, respectively. We choose 
the parallelogram of the fundamental periods as in Figure|3l As a base point Xo(fc) G 
Ffc we choose Xo{k) = (D]^{to{k)) where to{k) = T{k)/4. Let us notice that from (|24|)- 
(l26t it follows that for C > 1 we have 

Mi{to{k),k) = 1, S2{to{k),k) = -k', S3{to{k),k) = k. (80) 

Now, we consider four loops 

A{k), A'{k), y{k), y(fc):[0,l]^Ffc. 
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The loops A{k) and A'(fc) correspond to the real and imaginary periods, respectively 
(i.e. they correspond to the loops a and oc' in the parallelogram of the periods, see 
Figure |3l). 



T'ik)- 



to{k) 



h{k) 



sl(k) 



sKk) 



sl{k) 



h{k) to{k)+T{k) 



Figure 3: Parallelogram of periods with chosen paths. The points marked by ti(k), 
with i = 0, 1, 2, 3 are crossing points of the loops a, a', 13 and 13'. 

The loops Y{k) and y'(fc) corresponds to the 'shifted' real and imaginary periods, 
i.e. the loops /3 and (3' in the parallelogram of the periods. 

Remark 18. Above we use informal language. The correspondence between loops 
on r{k) and paths on the complex time plane can be viewed as follows. The map 

C3t — > 0{t,k) G r{k), 

is a covering map. For a loop a{k) on r{k) we obtain a path &{k) on C which is a 
lifting of cr{k) with respect to 0{-,k), i.e. a{k) is defined as such curve for which 

a{k) = 0{-,k)oa{k). 

These four loops cross at four common points Xi{k) = (D{ti{k),k), I = 0,1,2,3 
where ti{k) = to{k) + T{k)/2, t2{k) = to{k) + T'{k)/2 and = to{k) + T{k)/2 + 
T'{k)/2. Moreover, we have 

Miihik)) = 1, S2{h{k),k) = -k', S^{h{k),k) = -k, 
Mi{t2{k)) = 1, S2{t2{k),k) = +k', S3{t2{k),k) = +k, (81) 
Mi{t3{k)) = 1, S2{t3{k),k) = +k', S3{t3{k),k) = -k. 
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Thus, as k tends to 1, the points Xi{k) tend to u± and the loops A{k) and y(fc) ap- 
proach CI. We show that the loops A'{k) and y{k) tend to u±. In fact, for t = 
to{k) + It, T G M (i.e. along the loop ex.') from formulae (|24)l - ((26)l we deduce 
that S2{t,k) and S3(f, fc) are real while Mi(f, A:) — 1 is purely imaginary. If we put 
Ml - 1 = iMi in (|30l) we obtain 



2 



^M? - ^ + ^(S2 + CS2) = sl + Sl = 1, 



1 

2^ 



and thus 

It follows that for k = 1 



-Mi + cv^{l- Si) = a?k^. 



Ml + a?Sl = 0, 



but Ml and S2 are real, so we have Mi — 1 = and S2 = 0, i.e. the loop A!{k) tends 
to M+ as A; tends to 1. Similarly we show that A' {k) tends to W- as A; tends to 1. Four 



r(k)- 



rl{k) 



nik) 



Qfl 



0(2 



h{k)+T{k) 



Figure 4: Loop (54 is homotopic with loop a ■ a[ ■ (32^ ■ I5[ ^ ■ a-^^ , where (x = (X\- ai. 

points Xiik), I = 0,1,2,3 divide four loops A(k), A'(k), Y(k) and y'ik) into eight 
semi-loops Aj{k), Aj(/c), Yi{k) and y'lik) which correspond to eight semi-loops aj, a\, 
I3i and IS'^, I = 1,2 in the parallelogram of the periods. Of course we have A{k) = 
Ai{k) ■ A2{k), a = ai- 0.2, etc. We show that the fundamental group n{r{k), XQ{k)) is 
generated by closed loops which are appropriate compositions of these eight semi- 
loops. The fundamental group 7t{r{k),xo{k)) is generated by homotopic classes of 
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six loops with the base point at XQ{k): A{k), A'{k), and four loops cr/(fc) encircling the 
singular points S/ (k), I = 1,2, 3, 4. They satisfy the following condition: 

0-1 (it) • o-2(fc) • a3(k) ■ a^{k) = A(k) ■ A'{k) ■ A{k)-^ ■ A'{k)-^. 

We show that the loop cr/ [k) has the same homotopic class as an appropriate com- 
position of the semi-loops Ai{k), Yi{k) and y'lik). For example: 

[cr,{k)] = [A{k) . A[{k) . Y2{k)-' • y'lik)-' ■ M{k)-']. 

Let ^4 be the loop encircling T4(A:) and let 0-4 (fc) correspond to 64^. Then, we easily 
deduce that ^4 has the same homotopic class as a - a[- ^2 ^ " " '^1^' see Figure|ll 
Thus, we show that all generators of the fundamental group 7t{r{k),xo{k)) ap- 
proach n as A: tends to 1. □ 

Now, we are ready to prove our main result. 

Theorem 3. J/C ^ 1, A = B = 1, L = (0,0, 1), and 2£, ^ 3(1 - C), then system © 
considered on Ai^ does not admit an additional real meromorphic first integral functionally 
independent together with H and H5 in a neighbourhood of the phase curve Vi. 

Proof Let us assume that such meromorphic integral exists in a real neighbourhood 
of the phase curve Pi. Then we can extend it to a complex meromorphic first integral 
in a complex neighbourhood U oiVi. Then, by Lemma |8l we find Vj^ with k close to 1 
such that its fundamental group is generated by loops which lie entirely in U. From 
the Ziglin LemmaHJit follows that the monodromy group of the NVE ()4T|l possesses 
an invariant. But the NVE (|4T|l is a Fuchsian equation and thus if its monodromy 
group possesses an invariant, then its differential Galois group also possesses an 
invariant, see Theorem 3.17 in |0]- However, by Lemma|5lwe show that the identity 
component of the differential Galois group of (|47|), and thus the identity component 
of the differential Galois group of (|4T|l . is SL(2, C). It follows that it does not possess 
an invariant, see Example 2.11(b) from L4J. A contradiction finishes the proof. □ 

6 Comments and Remarks 

Although, as it is commonly believed, most systems are not integrable and inte- 
grable systems are extremely rare, the example considered in this paper shows that 
to prove the non-integrability one has to use rather involved techniques. Never- 
theless, a proof of non-integrability of a system gives, in some sense, a negative 
result — the true aim is to find a non-trivial integrable system. From this point of 
view, the reader can wonder why we did not investigate carefully the case of the 
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parameter values 2^ = 3(1 — C) for which the necessary conditions for integrability 
are satisfied. As a matter of fact, for some time we believed that for these parameter 
values the system is integrable. With the help of the computer algebra we tried to 
find a polynomial or rational first integral of the system but we failed. For fixed val- 
ues of C we numerically generated the Poincare cross sections of the system which 
evidently showed that the system is not integrable. Thus, our conjecture is that the 
system also is non-integrable for the case 2<^ = 3(1 — C). An analytic proof of this 
fact needs a separate investigation. 

For ^ = Theorem |3l tells us that the problem of a symmetric rigid satellite in a 
circular orbit is not integrable for all values of C G (0, 2) except C = 1. This problem 
was also investigated in [20; 21; 3] where a proof of the same fact is given. However, 
in all these references as a particular solution a heteroclinic orbit was chosen and in- 
stead of transformation (|43|) another one was used. This leads to a more complicated 
form of the NVE. 

7 Acknowledgements 

We would like to thank Delphine Boucher, Juan J. Morales-Ruiz, Jacques- Arthur 
Weil, Carles Simo, Michael F. Singer and Felix Ulmer for discussions and help which 
allowed us to understand many topics related to this work. As usual, we thank 
Zbroja not only for her linguistic help. For the second author this research has been 
supported by a Maria Curie Fellowship of the European Community programme 
Human Potential under contract number HPMF-CT-2002-02031. 

Appendix 

Let us consider a linear second order differential equation with rational coefficients 

w" + p{z)w' + q{z)w = 0, p{z),q{z) eC{z). (82) 

A point z = c G C is a singular point of this equation if it is a pole of p(z) or q{z). 
A singular point is a regular singular point if at this point p{z) = (z — c)p(z) and 
q{z) = (z — c)^q{z) are holomorphic. An exponent of equation ((82)l at point z = zq is 
a solution of the indicial equation 

p(p-l) +poP + '?0 = 0, po = P(c), qQ = q{c). 
After changing the dependent variable z ^ 1 /z equation (l82t reads 

w" + P{z)w' + Q{z)w = 0, P(z) = Ip (-\ + - Q(z) = 1^ (-\ . (83) 
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We say that the point z = oo is a singular point for equation (l82t if z = is a 
singular point of equation (l82b . Equation (l82b is called Fuchsian if all its singular 
points (including infinity) are regular, see |36: il5|l 

If one (non-zero) solution Wi of equation ((82)l is Liouvillian, then all its solutions 
are Liouvillian. In fact, the second solution IV2, linearly independent from ivi, is 
given by 

W2 = ivi I — exp 



w 



1 



Putting 



w = y exp 



into equation (l82t we obtain its reduced form 



y" = r{z)y, r{z) 



-q{z) + ^p'{z) + ^p{zf. 



(84) 



(85) 



This change of variable does not affect the Liouvillian nature of the solutions. For 
equation dSSt its differential Galois group Q is an algebraic subgroup of SL(2, C). 
The following lemma describes all possible types of Q and relates these types to 
forms of solution of (|85|l, see IitLEsIi . 

Lemma 9. Let Q he the differential Galois group of equation (|85)| . Then one of four cases 
can occur. 

1. Q is conjugated with a subgroup of the triangular group 



T 



a b 
fl-i 



aeC^bey, 



in this case equation ((85)l has an exponential solution, 
2. Q is conjugated with a subgroup of 



c 




U 



c 




in this case equation (|85)| has a solution of the form y = exp / w, where cu is algebraic 
over C(z) of degree 2, 

3. Q is primitive and finite; in this case all solutions of equation ((85)l are algebraic, 

4. Q = SL(2, C) and equation ((85)l has no Liouvillian solution. 
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When the first case occurs we say that equation dSSt is reducible. 

The Kovacic algorithm jl^ allows to decide if an equation of the form dSSt pos- 
sesses a Liouvillian solution. Applying it we also obtain information about the dif- 
ferential Galois group of this equation. Now, beside the original formulation of this 
algorithm ^ we have its several versions and improvements and extensions to higher 
order equations [U;JA;M;22;2Q;^M 

Here we present a part of the Kovacic algorithm which allows to decide whether 
(|85)l possesses a solution of the form y = exp / w, where cu is algebraic over C(z) of 
degree 1 or 2, or, in other words it gives an answer whether for equation ()85)l case 
1 or 2 in Lemma |9l can occur. We used this part of the algorithm in Lemma HI and 
LemmalZl As our NVE is Fuchsian, we present the algorithm adopted for a Fuchsian 
equation because it is simpler than for a general case. 

We write r(z) G C(z) in the form 

r(z) = s(z), t{z) G C[z], 

where s(z) and t{z) are relatively prime polynomials and t{z) is monic. The roots of 
t{z) are poles of r(z). We denote I' := {c G C | t{c) = 0} and I := I' U {oo}. The 
order ord(c) of c G I' is equal to the multiplicity of c as a root of t{z), the order of 
infinity is defined by 

ord(oo) := max(0,4 + degs — deg t). 

As we assumed, equation ((85)l is Fuchsian, so we have ord(c) < 2 of c G 1. For each 
c G 1' we have the following expansion 

r{z) = ^^ + o(-^), 
(z - c)^ \z- c J 

and we define Ac = a/I + 4flc- For infinity we have 

and we define Aoo = \/l + 4floo- 

Now we describe the Kovacic algorithm for the two cases mentioned. 

Case I 

^On the web page http://members.bellatlajitic.net/ jkovacic/lectures .html the reader 
will find lecture notes of J.J. Kovacic which contain an extended description of the algorithms with 
many remarks and comments concerning recent works on the subject. 
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Step I. For each c G Z' such that ord c = 1 we define Ec = {1}; if ord c = 2 

Ec := |^(l + A,),i(l-A,)|. 
If ord (oo) < 2 we put Eqo = {0, 1}; if ord (oo) = 2 we define 

Eoo := |^(l + Aoo),^(l-Aoo: 

Step II. For each element e in the Cartesian product 

E •= Eoo X 1^ Ec, 

C6l' 

we compute 

We select those elements e G E for which is a non-negative integer. If there 
are no such elements equation (|85)) does not have an exponential solution and the 
algorithm stops here. 

Step III. For each element e G E such that = n G Nq we define 

^(z) = E 



^r/ Z — C 



and we search for a monic polynomial P = P{z) of degree n satisfying the following 
equation 

P" + lcu{z)P' + {w'{z) + w{zf - r{z))P = 0. 



If such polynomial exists, then equation (l85t possesses an exponential solution of 
the form y = P exp / cu, if not, equation dSSb does not have an exponential solution. 

Case II 

Step I. For c G I' such that ord c = 1 we define Ec = {4}; if ord c = 2 

E, :={2,2(l + A,),2(l-A,)}nZ. 
If ord (oo) < 2 we put Eoo = {0,2,4}; if ord (oo) = 2 we define 

Eoo := {2, 2(1 + Aoo),2(l - Aoo)} n Z. 
Step II. If the Cartesian product 

E Eoo ^ E^;, 

C6l' 
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is empty then case 2 cannot occur and algorithm stops here. If it is not, then for 
e G £ we compute 

die) := 2 - - £ ec. 

We select those elements e G £ for which d{e) is a non-negative integer. If there are 
no such elements case 2 cannot occur and algorithm stops here. 
Step III. For each element e G £ such that d{e) = n ^ Nq we define 



6 = e{z) = 




and we search for a monic polynomial P = P{z) of degree n satisfying the following 
equation 

P'" + 39P" + (302 + 39' - 4r)P' + {9" + 399' + 9^- 4r9 - 2r')P = 0. 

If such a polynomial exists then equation dSSb possesses a solution of the form y = 
exp / cu, where 

11 P' 
cv^ - -ipcu + -il^' + -il>^ - r = 0, T^ = 9+—. 

If we do not find such polynomial, then case 2 in Lemma |9l cannot occur. 
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